Fermion-antifermion mixing in gravitational fields 
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Mixing of fermion and antifermion states occurs in gravitational interactions, leading to non- 
conservation of fermion number above temperatures determined by the particle masses. We study 
the evolution of a / , / system and calculate the cross sections for the reactions / ^^ /. Their values 
are identical in both directions. However, if / changes quickly into a lighter antiparticle, then the 
reaction symmetry is broken, resulting m an increased production of matter over antimatter. 

. ' Of the fermions that exist as free particles, baryons appear to satisfy, with good accuracy, a number conservation 

f~*) , law in the strong, electromagnetic and weak interactions. The same holds true, separately, for the numbers of the 

Cn| ■ lepton families. It is however thought that baryon and lepton numbers are not exactly conserved quantities. This is 

^ , sometimes invoked as a possible way to explain the value ns/n-y ~ 10~^ for the baryon to photon ratio observed in 

Q-i' the universe and which is at the heart of the baryosynthesis problem [l| . 

■^^ , It is shown below that fermion numbers need not be conserved in gravitational interactions. 

j«^ ' We consider external gravitational fields for which solutions of the covariant Dirac equation [2|-|4| 

^; [z7^(x)P^-m]*(a;)=0, (1) 

CJ ' and of other covariant wave equations [5|-l7| can be constructed. These solutions are exact to first order in the metric 
^ ] deviation 7^1, — g^y ~ rj^^, where r^^v is the Minkowski metric and can be calculated explicitly by means of path 
^ ■ integrals if the solution of the free wave equation (corresponding to the same particle of mass m) and ^^i^{x) are 
^0^ known. The notations are those of [J, Q, in particular V^ = V^ + ir^(a;), V^ is the covariant derivative and r^(a;) 
the spin connection. A semicolon and a comma are frequently used to indicate covariant and partial derivatives 
respectively. The first order solutions of ([1]) are of the form 
> 
O ■ *(a^) = T{x)^/j{x) , (2) 

fvq ' where i'{x) satisfies the usual flat spacetime Dirac equation 
^' 
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(lY'd^ - to) ipix) = , (3) 

and is represented here by the positive (particle) and negative (antiparticle) energy solutions '0(2;) — u(fc)e~*'^''^ and 
4>^^'{x) = w(fc)e*'^f^ respectively, where the spin-up ('f) and spin-down (j,) components of the spinors u and v obey 
the well-known relations [9| 

H = T^Vf , vi = j^u^ . (4) 



jrt The matrices 7'^ in ([S]) are the usual constant Dirac matrices. The operator T is given by 



f=--^(-*7'^(x)2?^-TO)e-**-, (5) 



where 



and 
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$T = ^s + $r 



^s{x)=f dz^Txiz), (6) 



^g{x) = -\f dz^ [7aA,M^) - IpxAz)] [(x" - z") fc^ - {x^ - 7J) fc"] -f ^ y" dz^7oA(^)fc" . (7) 

In ([6]) and (O , the path integrals are taken along the classical world line of the particle starting from a point P and 
$G ^ 0["ff,^) is the covariant Berry phase [l3|. It follows from (g]), 7^^ = 170717273 , {75,7'^} = 0, a^f^ = §[7", 7^] 
and 

r{x)^e%{x)r. r^(2;)=-|^"'^e:^e,3;M' [7^,^=0, (8) 



that if ^{x) = e-''""'''' fu is a solution of ([J), then *(i)(a;) = e^'^f^^fiz; also is a solution of (P), and fi = 75f75. It 
is useful to further isolate the gravitational contribution in the vierbein fields by writing e^ ~ (5^ + ft.^ . We therefore 
obtain 

f^fo + fG^^[{l- i^g) {m + 7"fca) -i{m + 7"fc,) $5 + (fc^/i^ + *G,a) 7"} , (9) 

where Tq = ^ (m + 7"A;q) and Tq contains the (first order) gravitational corrections. Similarly, we find 

Ti = Tio + ^iG = ^ {(1 - **g) (m - 7"fc„) - i (m - 7"^,) $5 - (A:/;/^^ + $G,a) 7"} , (10) 

where Tig = -^ (m — j^ka). Notice that i){x) and ^p'^^^x) as well as ^(x) and ^''^^•'(a;) describe particles of the same 
mass m as required by CPT-invariance. 

In general the eigenstates U^ = l/-\/2(u ± v) of 7^ and the eigenstates u and w of Tq are not the same and T, Ti 
mix u and u, thus producing // transitions. A similar mixing phenomenon occurs in the K'^ , K^ system where the 
CP eigenstates are mixed by the weak interaction Hamiltonian resulting in CP violation. It is the fermion number 
conservation that is violated in our work and the mixing is carried out by the gravitational contributions Tq and Tiq. 
This can be seen as follows. 

The time evolution of a /, / system in a gravitational field can be written in the form 

Mt)) = ao|*(t)) + M^^'^it)) = aof{tmt)) + /?oTi(t)|^«(i)) ^ «(i)|^(t)) + mH^'Ht)) , (H) 

where \ao\'^ + |/?oP — li from which we obtain 

a(t) = (V'|$(i))=ao(V'|T|'/')+/3o(^|Ti|V'(i));/?(i) = (V/''|$(i))=ao(^''^|T|^)+/3o(^'')|Ti|i/^(i)). (12) 



If at t = the gravitational field is not present, then Tq = , Tic = and ao = a(0) , /3o = /5(0). It follows from ^ 
that as the / , / system propagates in a gravitational field, transitions / ^ / can take place. 

In order to obtain the transition probabilities |a(i)p , |/3(t)P from (fT2|) in a concrete case, we choose for simplicity 



E + m f fn. \ __,t 



where f^ is the positive helicity eigenvector. The normalizations are (V'lV') = 1: where ('0| = (V'^|7°, (V'^^'^^IV' ) = ^1 
and (i/'IV''^'') — ('0'^''l^) — 0- In addition, we need explicit expressions of the metric components for the purpose of 
calculating T and Ti. We choose the Lense-Thirring metric 11|, represented, in its post-Newtonian form, by 

2 

700 = 2^ , 7y = 2(t>5ij , 7oi = ft-i = — (J A r)^ , (14) 



where 

GM . 4GA'fi?2 



-(2/,-x,0), (15) 



r hr-^ 

and M, R, uj ~ (0, 0, lu) and J are mass, radius, angular velocity and angular momentum of the source. This metric is 
non trivial, has no Newtonian counterpart and describes the physically significant case of a rotating source. By using 
the freedom allowed by local Lorentz transformations, the vierbein field to 0(7^^^) is 

e? = 0, e9 = l-0, ei^h,, e^ = (1 + 0) ^ . (16) 

Without loss of generality, we consider particles starting from z — —00, and propagating, with impact parameter 
b > R, along x = b ,y ^ in the field of the rotating source. We also set k^ = k and k^ = E. 

We can now return to (fT2)) . If originally the system is an antifermion, then ag = 0,/3o = 1, |'&(i)) = ri|-0'^^^) and 
from dTT]) we obtain 



a(i) = (V|fi|VW) 



2m 



,0„,0 , / 7,1,3 ' ^ 



Ehl^'> + ( -khl + ( — + fc ) <^(z) ) 7^ 



IV'^'^) , (17) 



where qi = k'^^ + ki, k'^ refers to {ip\ and go = 2i? because the external field is time-independent. The probability of the 
transition ip ~^ ip^ which, because of the coupling of spin to gravity, violates / number conservation, follows from 
(HH) and is 



V)<i 



Ht)? 
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E\ 



4>''{z). 



(18) 



In the approximation fc ^ to we obtain from (fT8|) P^(i)^^ — (3/4)^^^, while P.mi)^^ — {m/2k)'^(p'^ for k ^ m. The 
latter result leads in the limit fc — > to a well known infrared divergence that will be discussed below. 

If we choose ag = 1, /3o = 0, then |/3(i)p represents the probability for the inverse process ^ — ^ -0'^^^. We find 



P, 



j/>-s-i/)(i' 



\m\' = K^j^'Mmi' = K^^^'Mi'T.j'm' = \{m\^^''^)\ 
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(19) 



According to (J18p and (1191) , the transitions proceed in both directions with the same probability, with no contribution 
from the source's rotation when k is in the z-direction, as specified. 

With the help of (IT51) we can also calculate the cross section for the / — ^ / process which corresponds to a change 
|AL| = 2 for leptons and \AB\ = 2 for baryons. We get [III 



da 

dh 



Ek 



kjSiEf - E)dkf = Mf,\,p^'^)\\ 



where use has been made of the relations kfdkf = EfdEf , fc? 



Ej 



E^ 



(20) 
fc^, integration over Ef has 



been performed and the index / indicates momentum and energy in the final state. There is, of course, no energy 
transfer to the external field, since the field is static. We also find from (TTSl) 



MTiH^'ni 



(1)\|2 - 



2to2 ^ k ' 



m\ 



(21) 



where 0(fc) = —{GM/TT)KQ{bk) is the Fourier transform of (j){z) = —GM/y/z^ + lP, the momentum transfer is 
k =^ \kf — k\ — 2ksm{9/2) and 9 is the angle that kf makes with the z-axis. Integrating over the angles, we obtain 
from (P0|) the total cross-section 
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On account of the exponential decay of the Bessel functions, ([22| can be approximated by 

2 / ^ , r \ 2 
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(22) 



(23) 



for all physical values of bk. In (|23p . /3 and 7 are the usual Lorentz factors and k = m/Sj. It also follows from (PT|) 
and T = 7^X17^ that the transition amplitude ('(/''^^■'|T|-0) leads to the same cross-section ([2^ . From < /3 < 1, one 
obtains T > m/ k, = Tc and also Tc ~ 1.2 • W^K for neutrinos of mass to, ~ leV and Tc ~ 1.1 • 10^^ A' for nucleons. 
Expressing /3 and 7 as temperatures, we can re-write (|23p in the form 



1 ,GM 
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(24) 



where T = T /Tc. The transitions discussed are forbidden for T < Tc, have vanishing value for T — >■ 00 and diverge 
for T ^ 1 as expected. This divergence, already met in connection with (fT8| , arises as a consequence of the implied 
requirement that particles interact without the emission of gravitons (or photons in the electromagnetic case) [9(). 
This requirement can not in fact be realized physically. Processes in which gravitons with energy less than the energy 
resolution Ae of the apparatus are in fact indistinguishable from those in which gravitons of energy less or equal to 
Ae are also emitted and reabsorbed by the particles. The natural cutoff for the graviton energy is therefore Ae. It is 
also shown in the literature that when the gravitational 13 1 and electromagnetic 9] fields are quantized, all infrared 
divergencies disappear. 



Summary and discussion. It is possible for baryons and leptons to achieve the required energies in astrophysical 
situations. If, for the sake of numerical comparisons, we take GM/b ~ 0.1, which applies in the vicinity of a black hole, 
and /3 ~ 0.95, then the cross section value is cr ^ 8 • 10^^'^cm^ for leV neutrinos, while we find a ^ 8 ■ 10~^^cm^ for 
baryons of mass m ~ IGeV. It therefore follows from (P^ and 1^^ that, for the same value of /3, the baryon number 
is more likely to be conserved {a smaller) than the neutrino number, i.e. the S-violating processes are suppressed 
relative to those violating L. The actual value of a depends on (GM/b)^. It may therefore be useful to re-examine 
attentively some of the processes that take place close to compact, massive objects. 

Gravity behaves similarly on the cosmological scale and does not conserve the fermion number at lower momenta 
and higher temperatures. The values of T^ given above indicate that B conservation occurred in the universe at the 
time of primordial nucleosynthesis, while L conservation is relatively recent, thus resulting in the coexistence of two 
fermion populations at temperatures Tc^l ^ T < Tc^b, with the leptons agitated by the / ^ / transitions, and the 
baryons dormant (neglecting, of course, the effect of any other concomitant interactions). 

The cross section for / ^ / reactions is the same in both directions, hence the B- and i-violating processes do 
not produce more fermions than antifermions. The baryosynthesis problem can not therefore be solved using the 
gravitational non-conservation of fermion number without additional assumptions regarding the initial distributions 
of matter and antimatter, or the breakdown of thermal equilibrium. 

There is, however, an important, additional possibility suggested by (1231) . If, in fact, / of mass m produced in the 
reaction / — > / changes rapidly into an antifermion of mass m' < m, then the cross section for the inverse /' — > /' 
process can no longer re-establish symmetry between matter and antimatter. This reaction cross section does in fact 
increase, relative to the original one, by a factor a' /a ~ {m/m')^['y'{(3' — l//3'^)]^/[7(/3 — 1//?^)]^ which can be large, 
depending on {ra/m'Y' and the reaction kinematics. The result is therefore an increase in the production of matter 
over antimatter. Any further analysis of the mechanism's relevance to baryon asymmetry depends on m and m' and 
their experimental signatures and is not pursued here. 



Scalar-tensor theories of gravitation based on early efforts by a number of authors [lj| , [15| , [16| have been recently 



considered for the purpose of obtaining a varying gravitational constant that could in principle enhance the strength 



of gravity [19| , [17[ , [18| . Would the mixing mechanism discussed above be affected by the introduction of a running 
gravitational constant and concomitant scalar field? Unfortunately, no useful, first order tensor-scalar solutions of ([T]) 
exist at present. Nonetheless, a qualitative indication of the role played by a scalar field and running gravitational 
constant can be obtained from ([T]) by means of the conformal transformation 5^1, = G'A(a;)?7^i, which separates the 
scalar from the tensor component of the theory. In this simple case, ([T]) can be solved [20| , the operators T and Ti 
given by ^ and (|10l) become 

^' = i {("^ + 7"A:„) + I (lnGA(x))_„7"} (25) 

^' ^ i {("^-7-A:.) - I (lnGA(x))_„7"} (26) 

and still confirm the presence of the symmetry breaking mechanism. However, in view of the limitations on the 
structure of A(a;) discussed in [203, more precise conclusions must await the detailed study of a meaningful physical 
model. 

The introduction of back reaction terms discussed in ^8,] , while increasing the cross-sections over appropriate time 
intervals, does not alter the symmetry of the original distribution. The back reaction would however increase the 
cross-section of the process /' — >■ /' discussed above, and of a reaction for which CP violation has been observed, 
as in the time evolution of a K^K^ system, if the latter process occurred in the neighborhood of a gravitational 
source of appropriate strength. The model of back reaction discussed in [8j is based in fact on the introduction of a 
dissipation term in the wave equation of the particle propagating in an external gravitational field. It is then shown 
that a gravitational perturbation can grow rapidly, as with a fluid heated from below in which a small disturbance in 
the wave equation grows rapidly as soon as convection starts. The covariant wave equation of a kaon in the external 
gravitational field approximation has the solution f{x) = e~'^^^^''^'tfQ[x) where fQ{x) satisfies the equation [611 

{jrd^,d,+m^)^^{x)=Q. (27) 

It follows [8| that the addition to ((27|) of a dissipation term —2m£,doipo where ^ = ^ij^^^Y and e is a dimen- 
sionless, arbitrary parameter < e < 1, transforms the solution of the covariant wave equation into (p{x) — 
q-^'^g{x) ^m^xo i^^(^^^ ^j^(j |.]-^g overlap of the kaon eigenstates Kl and Ks, which is a measure of GP violation, into 
{Kl\Ks)g = {KL\Ks)e'^"^^^" . Then the increase in GP violation could become large even if the growth of the 
exponential term were restricted by competing effects [Sj. 
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